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University of Michigan 

By a well-known result of Bayer and Diaconis, the maximum 
entropy model of the common riffle shuffle implies that the number 
of riffle shuffles necessary to mix a standard deck of 52 cards is either 
7 or 11—with the former number applying when the metric used to 
define mixing is the total variation distance and the latter when it 
is the separation distance. This and other related results assume all 
52 cards in the deck to be distinct and require all 52! permutations 
of the deck to be almost equally likely for the deck to be considered 
well mixed. In many instances, not all cards in the deck are distinct 
and only the sets of cards dealt out to players, and not the order 
in which they are dealt out to each player, needs to be random. We 
derive transition probabilities under riffle shuffles between decks with 
repeated cards to cover some instances of the type just described. We 
focus on decks with cards all of which are labeled either 1 or 2 and 
describe the consequences of having a symmetric starting deck of 
the form 1,..., 1,2,..., 2 or 1,2,..., 1,2. Finally, we consider mixing 
times for common card games. 


1. Introduction. The connection between examples and concepts in prob¬ 
ability theory is a particularly close one. That examples derived from the 
question “How many shuffles mix a deck of cards?” have featured promi¬ 
nently in the development of the convergence theory for Markov chains by 
Persi Diaconis and others can be seen in this light. This article deals with 
riffle shuffling, which is the most common way of shuffling cards. 

There are 2” ways to cut a deck of n cards into two packets and then riffle 
them together since a card that ends up in the ith position can be dropped 
by either the left hand or the right hand. The maximum entropy model 
assigns equal probability to all these 2”’ riffle shuffles. More generally, the 
maximum entropy model assigns equal probability to all aP' a-shuffles, with 
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an a-shuffle being a way to cnt a deck into a packets and then riffle them 
together. Several equivalent descriptions of the a-shuffle have been given by 
Bayer and Diaconis [2] . The a-shuffle with a = 2 is also described by Epstein 
[5] who calls it the amateur shuffle. 

We will refer to elements of the group Sn of permutations of {1,2,..., n} 
as shuffles, tt € Sn and 7r(i) = j, then by convention the shuffle vr sends 
the card in the fth position to the jth position. The number of descents of 
vr is defined as the number of positions 1 <i <n at which 7r{i) > 7r{i + 1). 
Bayer and Diaconis [2] proved that the probability that an a-shuffle results 
in a shuffle vr with d descents is given ^). The picture below 

shows a 3-shuffle of 6 cards. 


2 0 110 0 



0 0 0 1 1 2 


The bottom line indicates that the 0th, 1st and 2nd packets in the cut have 
3, 2 and 1 cards, respectively. The top line indicates that the 1st, 2nd, 3rd, 
4th, 5th and 6th cards in the shuffled deck are dropped from the 2nd, 0th, 
1st, 1st, 0th and 0th packets, respectively. If the numbers are ignored, the 
arrows alone depict a shuffle. In a depiction of a shuffle such as the one 
above, a descent corresponds to a crossing between arrows that originate at 
adjacent positions. The shuffle depicted above has 2 descents, and therefore, 
according to Bayer and Diaconis [2] , the probability that an a-shuffle results 
in the shuffle depicted above is ^ ■ 

In nearly all of the literature on card shuffling, the cards in a deck are as¬ 
sumed to be distinct. We allow cards to be indistinguishable. In our notation, 
both 1,1,2,1 and 1^,2,1 denote the deck with two cards labeled 1 above a 
card labeled 2 above a card labeled 1. Let ai,a 2 , ■ ■ ■ ,an be a deck. When 
it is shuffled using vr G Sn, the deck obtained is Ojr-i(i),a^-i( 2 ), • • • ,a^-i(„). 
We define 7r(Di; D 2 ) as the set of shuffles n G Sn such that vr applied to Di 
results in D 2 . The descent polynomial of the shuffles from the starting deck 
Di to the ending deck D 2 is defined as 

^des(7r) 

7rG7r(Di;D2) 

where des(7r) is the number of descents in vr. For example, the descent poly¬ 
nomial of shuffles from 1,1,2,2 to 1,2,2,1 is 2x -|- 2x^. 

If the descent polynomial of the shuffles from a deck Di with n cards to a 
deck D 2 is the probability that an a-shuffle of Di results in the 
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deck D 2 is 


( 1 . 1 ) 



a + n — d — 1 
n 


a formula obtained by summing over all the shuffles in tt{Di ; D 2 ). The system 
of linear equations (1.1) can be inverted to obtain 



( 1 . 2 ) 



for 0 < d < n. It is possible to pass back and forth between the transition 
probabilities pa and the descent polynomial of Tr{Di;D 2 ) using (1.1) and (1.2). 

In Section 2 we deduce efflcient recursions for the descent polynomial 
from the starting deck Di to the ending deck D 2 when either Di or D 2 
is a sorted deck of the form , 2”^,..., . We derive a formula for the 

transition probabilities when Di = 1,2,... ,k,x'^ without using the descent 
polynomial. Sections 3 and 4 consider the starting decks Di = 1^,2” and 
Di = (1,2)”'. Section 5 summarizes mixing times for card games obtained 
using results in the preceding sections and Monte Carlo simulations. 

Although decks with repeated cards do not seem to have been considered, 
the work of Diaconis, McGrath and Pitman [4], Fulman [7] and Lalley [11] 
on cycle decompositions, and of Fulman [8] on increasing subsequences are 
in a somewhat similar vein. The thesis of Reyes [13] has new results, as well 
as many references related to other types of shuffles. 

2. Transition probabilities. We begin with a recursive algorithm to ob¬ 
tain the descent polynomial of shuffles from 1”^, 2”^, • • •, h”'* to a deck D 

which has the same ni + n 2 H- \-nh cards but in a different order. Each of 

the numbers ni,n 2 , ■ ■ ■ ,nh is a positive integer. The transition probability 
under an a-shuffle can be obtained using the descent polynomial and (1.1). 

We assume the starting deck to be 1”^, 2”^,..., /i”'*, which is in sorted 
order. We denote by D{i,c) the position of the ith card labeled c in D. 
For example, \i D = 1, 2,1,1,2,2, then D{2, 1) = 3. The deck obtained from 
D by keeping the cards labeled 1,2,..., e and by discarding cards with all 
other labels will be denoted Hie. Similarly, the deck obtained from D by 
keeping the cards with labels and discarding other cards will be 

denoted by Dj^. We assume l<---<e<f<---<h and that there is 
no card whose label is in-between e and /, or, equivalently, / = e -|- 1. Let 

= ni -F n 2 -I-h Uh, Ni = ni + n 2 -\ -h Ug and A ^2 = Then N, 

Ni and N 2 equal the number of cards in D, Hie and Dj^, respectively. 
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Consider the set of all shuffles tt from the sorted deck , 2”^,..., /i”'* 
to D such that '7r(l) = D{i, 1) and 7r{N) = D{j,h), where 1 <i <ni and 
1 < j < The number of these shuffles with d descents is set equal to the 
coefficient of to define the polynomial pij{x). 

To obtain a recursion for consider the set of shuffles from the 

sorted deck , 2”^ ^^ e""= to Die and the set of shuffles from the sorted 
deck /"■/,..., /i”'* to Dfh- Define qi^k{x), for 1 < i < ni and 1 < k < Ug, as 
the polynomial in which the coefficient of x'^ equals the number of shuffles vr 
with d descents belonging to the first set which satisfy 7r(l) = Die{i, 1) and 
7r(A^i) = Dig{k,e). The polynomial rij{x), for 1 <l <nf and 1 < j < rik, is 
defined similarly, with the coefficient of equal to the number of shuffles tt 
with d descents in the second set which satisfy 7r(l) = Dfk{l, f) and 7r(A^2) = 
Dfh{j,h). Then the following recursive relationship holds: 

(2.1) pij{x) = ^qi^k{x)rij{x)x^^^'^'>. 

k,l 

The indices k and I vary over 1 < k < rig and 1 <l < rif. The exponent £{k, 1) 
is 0 if D{k, e) < D{1, f) and 1 if D{k, e) > D{1, /). 

To prove (2.1), we consider a bijection between 7r(l"'i, 2”^,...,/i”'*; D) 
and 7r(l”i,..., e”®; Die) x 7r(/"^,..., K^^^Djh). Let the shuffle vr map to the 
pair of shuffles (vri,7r2) under this yet to be defined bijection. If position 
i is occupied by a card labeled 6 in the starting deck , 2”^,..., /i”'* and 
7r(i) = D{j,6), then 7ri(i) = Dig{j,6) if 1 < <5 < e and 7r2(i - iVi) = Dfh{j,S) 
if / < 5 < /i) by definition of the bijection. To complete the proof of (2.1), 
we relate the number of descents of vr to the number of descents of vri and 
vr 2 . The number of descents of vr equals the sum of the number of descents 
of vTi and 7r2 if vr(A^i) = D{k, e), 7r(A^i + 1) = D{1, f) and D{k, e) < D{1, /). 
However, if D[k,e) > D{l,f), the sum must be incremented by 1. 

The base case of the recurrence (2.1) occurs when the starting deck has 
cards of only one type. Take this deck to be 1”. The coefficient of x'^ in Pij{x) 
is then equal to the number of shuffles tt G Sn with vr(l) = i, vr(n) = j and 
with d descents. The number of shuffles it G Sn with d descents is dehned 
as the Eulerian number (^) [9]. Given a permutation of {l,2,...,n — 1} 
with d or d — 1 descents, the number n can be inserted in d + 1 or n — d 
places, respectively, to obtain a permutation of {1,2,..., n} with d descents. 
Thus, as shown in [9], consideration of the insertion of the number n into a 
permutation of the numbers 1,2,..., n — 1 gives the recurrence 

(d)=‘' (!!)=“ 
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The modified Eulerian number (”)j is defined as the number of vr G 5^ with 

7r(l) = i and d descents. If d = 0, (^ )i is 0 if z > 1 and 1 if i = 1. Consideration 
of the insertion of n into a permutation of the numbers 1,2,..., n — 1 that 
begins with i gives the recurrence 


(2.3) 



(d+ 1) 


n — 1 
d 




n — 1 
d-1 


if n = i, d > 0. 


if re > z, d > 0, 


If re = z = 1, (^ )i is equal to 1 if d = 0 and equal to 0 if d > 0. The modified 
Eulerian number {^)ij is defined as the number of vr G 5^ with vr(l) = z, 

T^in) = j and d descents. If d = 0, (^ )ij is 1 if z = 1 and j = n but 0 otherwise. 
Eor d > 0, the following recurrence can be derived: 




n> i,n> j, 


(2.4) 


re — 1 
d- 1 


n-j 



if re = z, re > j, 
if re > z, re = j. 


If re = z = J = 1 and d > 0, = 0. Using (2.2), (2.3) and (2.4), the poly¬ 

nomials pij{x) can be formed in the base case. 

The descent polynomial of shuffles from , 2”’2,..., d”'* to D is obtained 
as the sum of the polynomials pij{x) over 1 < z < rei and 1 < j < Uh- 

We now turn to the descent polynomial of shuffles vr from D to the sorted 
deck , 2”2,..., . We first consider the occurrence of descents between 

7r(/c) and 7r(A: -|-1) when the positions k and k + 1 are occupied in D by cards 
with different labels. There will be a descent if and only if the label of the 
card at k is greater than the label of the card at k + 1. Thus, the number 
of descents of this type is the same for every shuffle from D to the sorted 
deck and is equal to the number of places where a card with a greater label 
immediately precedes a card with a lesser label in the deck D. This quantity, 
which may be denoted by des(Z)), is called the number of descents in D and 
is extensively studied in [10] and [12]. 

We next consider descents between TT{k) and 7r(A: + l) only if both positions 
k and k + 1 are occupied in D by cards with the label c. The cards at k and 
k + 1 both have the label c if and only if A: = D{i, c) and k + 1 = D[i + l,c) for 
some integer i, 1 < z < rzc. To facilitate the counting of this type of descent, 
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denote the generating polynomial of the Eulerian numbers by 

r]n{x) [10]. If we pay attention only to cards with label c in the deck it will 
look like * * ccc* **cc**c with blocks of c’s separated by cards with different 
labels. Assume that the lengths of these blocks are given by mi, m 2 , ■ ■ ■, m^, 
with 7 being the number of blocks. Then ric = mi + m 2 + ■ • ■ + m^. Let 

C 

jc = 

and let ic = jc — n^c + 1- If tt is a shuffle from D to the sorted deck, then 
ic<T^[D{i,c)) < jc must hold for 1 < i < ric. The ric integers in [ic,jc] can be 
divided into sets of mi,m 2 , ■ ■ ■ ,m^ in ricl/{mi\m 2 l ■ ■ ■ m^\) ways. For each 
such division of these ric integers into sets, there are mi!m 2 ! • • • m-f\ ways of 
assigning values to 7r{D{i,c)), 1 <i< ric, such that a number assigned to a 
position within the hrst block of cs is in the first set and so on. The coefficient 
of x'^ of the polynomial f?mi( 2 ;)hm 2 (^)' ’ is equal to the number of 

these assignments which have d descents. Therefore, the coefficient of of 
the polynomial 

71 ^ 

(2.5) Pc{x) = -j- ^{x)Vm 2 (x)--- (x) 

is equal to the number of assignments with d descents out of the ricl assign¬ 
ments of integers in [ic,jc] to TT{D{i,c)), 1 < i < ric- As intended, (2.5) counts 
the descent between 7r(fc) and 7r{k + 1) if and only if cards at positions k 
and k + 1 in D both have the label c. 

To hnd the descent polynomial of shuffles from D to the sorted deck, 
note that the occurrence of a descent between n^k) and 7r{k + 1 ), with cards 
labeled c at positions k and /c -|- 1 in T), is completely independent of the 
occurrence of a descent between 7r(l) and 7r(l + 1), with cards labeled d at 
positions I and 1 -|- 1, if c 7 ^ d. Moreover, there are always des(D) descents 
in a shuffle tt from D to the sorted deck that correspond to positions k and 
k + 1 occupied in D by cards with different labels. Therefore, the descent 
polynomial is given by 

( 2 . 6 ) x'^‘^^‘^^'^pi{x)p2{x)---ph{x), 

where the Pi{x) are dehned by (2.5). 

If the deck D is any permutation of the multiset {l"'i, 2"'^,...,/i”'*}, 
( 2 . 1 ) and ( 2 . 6 ) make it possible to hnd the descent polynomials of shuf- 
hes from the sorted deck to D and of shuffles from D to the sorted deck in 
polynomial time. The descent polynomial of shuffles between decks neither 
of which is sorted will be considered in later work. 

In the rest of this section, we turn to theorems about transition proba¬ 
bilities between decks under an a-shuffle which do not use the descent poly¬ 
nomial. Let ai,a 2 , ■ ■ ■ ,an be one of the a"' integer sequences with 0 < a* < a 
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for 1 <i <n. This sequence can be sorted to < 0*3 < • • • < in a stable 
manner and the permutation ii,i 2 , ■ ■ ■ ,in of { 1 , 2 ,...,re} is uniquely deter¬ 
mined since we require ij < ij+i if . Associate the shuffle tt € Sn 

with 7r(A;) = for 1 < A; < re with the sequence oi, 02 ,..., a^. Then the uni¬ 
form distribution on the sequences induces the a-shuffle distribution on 
Sn [2]. This description of the a-shuffle is used in Theorems 2.1 and 2.2. 

Theorem 2.1. Among all decks D that are permutations of the multi¬ 
set {l”i, 2 ” 2 ,..., h^^}, the transition probability under an a-shuffle from the 
sorted deck , 2"^ ^^ h'^^ to D is greatest for D = , 2"^ ^^ h'^'^ and 

least for D = hJ^^^,... ,2'^-^ 

Proof. Assume the sorted deck to be 1”, 2”'. The proof for more general 
sorted decks is similar. 

Let oi, 02 ,..., a 2 n be a sequence with 0 < a* < o for 1 < i < 2re. If D{i, 1) = j, 
define a* = aj, and if D{i, 2) = j, define Pi = aj. For example, if H = 1,2,1,2,1,2, 
then 

0 . 1 , 02 , as, 04, 05 , qq = ai, Pi,a 2 , P 2 ,ois, Ps- 

For the sequence to induce a shuffle from 1"',2” to D, each Oi must be less 
than or equal to each /3j. In addition, each a* must be strictly less than all 
the P's that precede it in the sequence. For example, if H = 1,2,1,2,1,2, the 
inequalities are 

max(Q;i,Q!2,a3) < min(/3i,/32,/33), 0:2 < Pi, as<Pi, as < P 2 . 

If H = 1”,2”', it is enough if each a is less than or equal to each p. If 
H = 2”, 1”, each a must be strictly less than each /3. Therefore, the number of 
sequences that induce a shuffle from 1"', 2"' to D is greatest for H = 1"', 2” and 
least for D = 2^,V'‘. The statement about transition probabilities follows. 

□ 


Theorem 2.2 below generalizes Theorem 3 of [2] and their proofs use 
similar arguments. Similar arguments can also be found in [ 6 ] and [10]. 


Theorem 2.2. Let the deck D he a permutation of the multiset {1,2,..., 
/i,x”}. Let the number of cards labeled c, 1 <c<h, that are not preceded by 
a card labeled c — 1 in D be equal to r. Let the number of cards labeled x 
that precede the card labeled h in D be equal to 1. Then the probability that 
an a-shuffle applied to the sorted deck 1,2,..., /i, x"' results in D is 


1 


a—1 




m=r—l 


f m — r -\- h\ , .w. .n-l 

j{a-m-l) (a-m) , 


where if I = 0 and m = a — 1, {a — m — lY must be taken to be I. 
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Proof. Let ai,a 2 ,..., au+n be an integer sequence with 0 < a* < a for 
1 <i <h + n. If D{1, c) =i, 1 <c<h, define Uc = fli- If D{j, x) = k, 1 <j < 
n, define Pj = Ofc. For the sequence ai,a2, ■ ■ ■ ,ah+n to induce an a-shuffle 
from the sorted deck to D, we require 

(2.7) ai<a 2 <-- -<ah< min(/3i,.. .,/?„)• 

In addition, the inequality ac-i < etc, 2 < c<h, must be strict if the card 
labeled c in L) is not preceded by the card labeled c — 1. Therefore, exactly 
r — I inequalities between the as in (2.7) are strict. Further, at least I of the 
Pi, the ones with 1 <i < I, must be strictly greater than ah- 

The number of solutions to (2.7), with the additional conditions described 
below it, can be counted by allowing a/^ = m to vary from r — 1 to a — I. Given 
m, the number of ways to pick the Uc, 1 <c<h, can be counted as follows. 
Start with m “jumps.” Allocate r — 1 of these jumps to the inequalities in 
ai < a 2 < ■ ■ ■ < Oih-i P ni that must be strict. The remaining m — r + 1 
jumps can be assigned to h positions, namely, the position before a\ and 
the h — 1 inequalities, in ) ways. The value of a*, for 1 <i <h — 1, 

is equal to the number of jumps preceding it. The number of ways to pick 
the /3s is {a — m — l)^(a — The formula for the transition probability 

from the sorted deck to D follows. □ 

3. Starting deck !”■, 2”'. The probability distribution over decks that 
are permutations of the same multiset of cards under an a-shuffle can be 
obtained from ( 2 . 1 ) or ( 2 . 6 ) if either the starting deck or the ending deck is 
in sorted order. The total variation distance from the uniform distribution 
is a sum over all possible decks and its calculation can therefore involve a 
very large number of terms. However, the calculation becomes simpler if it 
is recognized that the transition probabilities are the same for whole classes 
of decks. In the case where all n cards have different labels, the transition 
probabilities depend only upon the number of descents in the shuffle and, 
hence, the n! decks fall into only n equivalence classes. In this section we 
investigate this type of equivalence relationship when the starting deck is 
1 ^, 2 ". 

In this section and the next, we use a, /3 and 7 to denote sequences of 
I’s and 2’s that stand for segments of a deck of cards. The number of entries 
in the sequence a is denoted by |a|. The sequence obtained by reversing the 
order of a and then replacing each 1 by 2 and each 2 by 1 is denoted a*. For 
example, if a = 1,2,2,2,1,1, then a* = 2,2,1,1,1,2. A total of (^^) decks 
can be obtained by rearranging the cards of 1”, 2"'. The equivalence relation 
R on that set of decks is defined as follows. The deck Di = aP'y is ii-related 
to D 2 = a/ 3*7 if |a| = I 7 I, and the number of I’s and the number of 2 ’s in 
P are equal. For example, 1,2,2,1 is ii-related to 2,1,1,2 and 1,1,2,2,1,2 
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is i2-related to 1,2,1,1,2,2. The equivalence relation is obtained by taking 
the transitive, reflexive closure. For example, the decks 1,2,1,2,2,1,2,1 and 
1 , 2 , 2 , 1 , 1 , 2 , 2,1 and 2 , 1 , 1 , 2 , 2 , 1 , 1,2 and 2 , 1 , 2 , 1 , 1 , 2 , 1,2 are all in the same 
equivalence class. 

Theorem 3.1. If Di is R-related to D2, the transition probability from 
I"", 2 ” to Di is equal to the transition probability from 1 ”, 2 "' to D 2 under an 
a-shuffle for any a. 

Proof. It is sufficient to consider Di = a^'y and D 2 = ajl*^ with |q!| = 
I 7 I and with equal number of I’s and 2’s in /3. It is enough to show that 
the descent polynomial of shuffles from I"", 2 "' to Di is equal to the descent 
polynomial of shuffles from 1”,2” to D 2 . We will construct a bijective map 
from 7 r(l”, 2”; Df) to 7 r(l"', 2"; D 2 ) such that a shuffle maps to another shuffle 
with exactly the same number of descents. 

The number of I’s in a is equal to the number of 2’s in 7 since the 
number of I’s and 2’s are equal in a/Jq and in [5. Similarly, the number of 
2’s in a is equal to the number of I’s in 7. Let pi < ■ ■ ■ <Pa and qi < ■ ■ ■ < qb 
and ri < • ■ • < Tc be the positions of Di that correspond to a and f3 and 
7, respectively, that are occupied by I’s. Similarly, let p'l < • ■ ■ < p'c and 
qi< • ■ ■ < q'b and < ■ • • < J'a be the positions of Di that correspond to a 
and (3 and 7, respectively, that are occupied by 2’s. 

Define f{x) = 2n — x + 1. The map / reflects a position in a deck of size 
2 n about its center; for example, the first position is reflected to the last 
position. The positions occupied by I’s in D 2 = q;/ 3*7 are 

<Pa< f{q'b) <■■■< f{q'i) < n < • • • < rc, 

where ps and rs correspond to a and 7 , and f{q)s correspond to (3*. In Di, 
the position q[ is occupied by a 2. When j3 is reversed that 2 is moved to 
the position f{q'j) and then it is replaced by 1 to form (3*. This explains the 
central block of fiq'fj's above. Similarly, the positions occupied by 2’s in D 2 
are 

p'i<---<Pc< fiqt) <■■■< f{qi) <r[<---<r'^, 

where the p's correspond to positions in a, f{q)'s to positions in (3* and r'’s 
to positions in 7 . Note that each p or p' is less than each q or q', which is 
less than each r or rh 

Let TT G S 2 n be a shuffle from 1”',2"' to Di. Then the numbers 

7 r(l), 7 r( 2 ),..., 7 r(n) 

must be an arrangement of the positions occupied by I’s in Di. Similarly, 
the numbers 


7r(n + l),7r(n + 2),... ,7r(2re) 
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must be an arrangement of the positions occupied by 2’s in Di. The map 
to a shuffle from to D2 is based on two cases. In the first case, we 

assume that not both 7 r(n) and 7 r(n + 1) correspond to positions in /3. The 
shuffle vr* from I”, 2 ” to D2 that vr maps to is dehned as 

7r*(i) = 0(7r(i)), 

for 1 < i < 2n, where (j){-) will now be defined. First, we define =Pi, 
(p{ri) = Tj, (t>{p'j) =p'i and = r[. In addition, we define 4> as 

Qi ^f{Qb) Q2^f{qb-i) ■■■ qb^fiq'i) 

q'l ^f{qb) q'2^ f{qb-i) ••• q'b^fiqi)- 

This definition maps the qs to the f{q')^s and the g'’s to the /((?)’s and, 
therefore, vr* is a shuffle from 1^,2” to D 2 . Further, x < y if and only if 
4>{x) < 4>{y), except when x is a q and y is a q' or when x is a q' and y is a 

q. However, in the arrangement 

7 r(I),... , 7 r(n), 7 r(n + I),... , 7 r( 2 n), 

a q and q' can occur in consecutive positions as 7r(n) and 7r(n+ 1), and in no 
other way, and we have assumed that not both of those positions correspond 
to p. Therefore, the above arrangement has the same number of descents as 

(/>(7r(l)),..., (/>(7r(n)), (/>(7r(n + 1)),..., (?i(7r(2n)) 

and vr* has the same number of descents as vr. 

The other case is when vr(n) is a g and vr(n + 1) is a q'. Then we dehne 

<p{qi) = f{qi), 4 >{q'i) = fiq'i), Pi ^ ^'a-i+i and We map vr to 

vr*, where vr* is dehned as 

vr*(f) = (j){7r{2n — v + 1 )), 

for 1 < i < 2n. It can be verihed that vr* is a shuffle from 1"',2” to D 2 . 
Further, x < y if and only if 4>{x) > 4>{y), except when x is a p and y is a p', 
or X is a p' and y is a p, or x is an r and y is an r', or x is an r' and y is an 

r. In the arrangement, 

7 r(l),..., vr(re), vr(n + I),..., vr(2re), 

a p and p' or an r and r' can occur in consecutive positions only at vr(n) and 
vr(n + I). However, we have assumed that 7 r(n) is a 5 and that 7 r(n + 1) is a 
q'. Therefore, every descent in the above arrangement becomes an ascent in 

0(7r(l)),..., (/>(vr(n)), (/>(vr(n + 1)),..., (?i(vr(2n)) 

and every ascent becomes a descent. This arrangement is reversed to dehne 
vr*(l),..., vr*( 2 n) which changes the ascents back into descents and, there¬ 
fore, the number of descents in vr* is equal to the number of descents in 
vr. 
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Finally, we need to show that the map defined above is a bijection. A 
shuffle from 1 ^, 2 ” to D 2 can be mapped to a shuffle from 1 ”, 2 ” to Di using 
the same procedure as above. The resulting map is the inverse of the above 
map because (po <f) is identity in both cases above. □ 

It is natural to ask if the equality of the descent polynomials of the shuffles 
from I”, 2”’ to Di and D 2 implies that Di is ii-related to 02 - We have 
checked that this is indeed so for n = 1,2,3,4,5, 6 ,7. The theorem below 
counts the total number of equivalence classes under the relation R. 

Theorem 3.2. The number of equivalence elasses under R is equal to 
the Catalan number Cn+i) ■ 

Proof. We describe a method to hnd a unique representative for each 
equivalence class and then count the number of unique representatives. The 
function f{x) = 2 n — x + 1 reflects positions with respect to the center of 
the deck as before. In this proof, we refer to f{x) as the reflection of the 
position X. A position and its reflection either both he in (3 or lie outside it, 
since |a| = I 7 I. Consider the positions x = re + l,n + 2 ,..., 2n. Assume that 
X and f{x) both lie inside [3. If the positions x,f{x) are occupied by 1,2 
in the deck D, reversal of j3 changes it to 2,1, and the replacements of I’s 
by 2’s and 2’s by I’s changes it back to 1,2. Similarly, application of the 
basic rule that generates the relation R does not change D in the positions 
X, f{x) if those positions are occupied by 2,1. However, if they are occupied 
by 1 , 1 , that becomes 2,2 when the rule is applied using a (3 large enough to 
include positions x and f{x). Similarly, 2,2 becomes 1,1. With each position 
X, we associate the symbol if positions x,f{x) are occupied by 1,2, the 
symbol ” if occupied by 2 , 1 , the symbol 1 if occupied by 1 , 1 , and the 
symbol 2 if occupied by 2,2. The deck as a whole is coded as the list of 
symbols associated with positions n + 1 through 2n. For example, the deck 
1 , 1 , 2 , 1 , 2 , 2 , 2 , 1 , 1 , 2 , 1 , 2 , 1,2 is coded as +, +, 2 , 1 , 2 , 1 , -. 

The +s and — s never change when the rule that generates the relation 
R is repeatedly applied with possibly many different choices of (3. They are 
ignored in much of the rest of this proof. We can find the /3’s which lead 
to a nontrivial application of the rule to generate the relation R using the 
code for the deck D as follows. We traverse the code from left to right, and 
record the excess of I’s over 2’s. For example, for the code —, —, 2,1,2,1, +, 
this excess is —1 after the first 2 is passed, then becomes 0 , and then — 1 , 
and then 0. The rule for generating R can be applied whenever this excess 
becomes 0. If the excess becomes zero, after traversing i symbols in the code, 
the corresponding f3 in the deck is a segment of 2 i cards extending from 
position n — i + 1 to position n + i. When the rule is applied, the I’s become 
2’s and the 2’s become I’s among the hrst i symbols of the code. If the excess 
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becomes 0 after i symbols and then again after i-\-j symbols, the application 
of the rule with a /3 of length equal to 2 i followed by another application 
using a /3 of length 2 {i + j) changes the code for the deck only between the 
(i + l)st and the jth symbol. Among these symbols, the I’s change to 2’s and 
the 2’s change to I’s. By applying the rule with judicious choices of /3, it is 
possible to obtain a single code in which the excess never becomes negative. 
For example, the code —, —, 2,1,2,1, + can be converted to —, —, 1,2,1,2, +. 
We use such codes as unique representatives of equivalence classes of decks. 

Assume that in such a code, there are k symbols equal to 1 and k symbols 
equal to 2. Then there must be n — 2A; symbols equal to + or —. If the + and 
— symbols are ignored, and each 1 is substituted by a ( and each 2 by a), we 
obtain a valid arrangement of parentheses of length 2k. The number of valid 
arrangements of parentheses of length 2k is well known to be the Catalan 
number For each assignment of I’s and 2’s to the 2k positions, 

the other positions can be filled with symbols + and — in 2”“^^ ways. The 
2 k positions that are assigned either the symbol 1 or the symbol 2 can be 
chosen in (^) ways. Therefore, the total number of equivalence classes is 

given by Yl,k \-+i Ck)i 2 k)- This sum can be simplified: 

^^ 2 k\ ( n \ ^^ns^^(n/ 2 \ ({n-l)/ 2 \ 

Z^k + l\ k ) \2k) ^k + 1 \ k jy k ) 

2- ^/n/2 + l\ /(n-l)/2\ 

re/2 + 1 k + 1 J \ k ) ' 

The first equality above uses (5.35) in [9]. The proof may be completed using 
the binomial identity Y.k (m+fc)(p+fc) = integers l,m,p and I > 0. 

The cases with re even and odd have to be considered separately. □ 

4. Starting deck (1,2)”. The equivalence relation R in this section is 
different from the one considered in the previous section. In this section 
Di = afj'y is i?-related to D 2 = a/3*7 if f3 has the same number of I’s as 
2’s. The additional condition |a| = lyj is no longer required. The decks 
Di are all permutations of {1”,2"'}. The equivalence relation is obtained 
by taking the transitive, reflexive closure. The equivalence class containing 
1,2,1,2,2,1,2,1 has five other decks in it. 

Theorem 4.1. If Di is R-related to D 2 , the transition probability from 
(1,2)"' to Di is equal to the transition probability from (1,2)" to D 2 under 
an a-shuffte for any a. 

Proof. It is sufficient to consider Di = ct/Sy and D 2 = q;/3*7 with equal 
number of I’s and 2’s in (5. We will construct a bijective map from '?r((l, 2)"; Hi) 
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to vr((l, 2)”; D2) such that a shuffle maps to another shuffle with exactly the 
same number of descents. 

Let TT e 7r((l, 2)"'; Di). Let 7r(2i — 1) = a*, 1 < i < n, and 7r(2i) = hi,l <i <n. 
The number of descents in vr is equal to the number of descents in the 
arrangement ai, 61,02,62 ,..., To facilitate the proof, we depict this 

arrangement in the following way: 


rti 02 


a 


n 



bx b-2 


b 


n 


In the deck Di each position a* is occupied by a 1 and each position 6* is 
occupied by a 2 [because tt is a shuffle from (1,2)” to Di]. We assume that 
/3 begins at the {i + l)st position and ends at the (i + 2j)th position. In the 
depiction above, circle all the a* and hi that do not correspond to f3, that 
is, circle an a* or a bt if it is less than i + 1 or greater than i + 2j. When tt 
is mapped, the circled numbers will stay fixed. The uncircled numbers form 
segments that run from a circled number to another circled number (or they 
might begin or end at ai or 6„). These segments of uncircled numbers are 
of four types according as they either begin or end at the top line or the 
bottom line: 





The third type of segment has one more uncircled position in the bottom 
line, corresponding to a position occupied by a 2 in Di, than in the top line. 
The fourth type of segment has an extra uncircled position in the top line, 
corresponding to a position occupied by a 1 in Di. Since the number of I’s 
in [5 is equal to the number of 2’s, the number of uncircled positions in the 
top line must be equal to the number of uncircled positions in the bottom 
line. Therefore, the number of uncircled segments of the third type must be 
equal to the number of uncircled segments of the fourth type. 

To map TT to a shuffle from (1,2)” to D 2 , we will modify the uncircled 
segments and insert them back in-between the circled numbers in the original 
arrangement of a* and 6*. We define a map / from the uncircled positions, 
that is, the positions that correspond to /3, back to the the uncircled positions 
as follows: 

i + 1 ^ i + 2j, i + 2^i + 2j — 1, ..., i + 2j — l—^i + 2, i + 2j ^ i + 1. 
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If i + 1 < X < i + 2j and the position x in Di is occupied by 1 (or 2), the 
position f{x) in D 2 will be occupied by 2 (or 1). If Op, bp, Op+i, 6p+i, ... ,bq is 
an uncircled segment of the first type, it will be modified to f{bq),f{aq),..., 
/(op+i), /(6p), /(up). In the deck D 2 , each position f{bi) [or /(a*)], p<i<q, 
is occupied by 1 (or by 2). Therefore, the modified segment is also of the 
first type. However, when an uncircled segment of the third (or fourth) 
type is modified in this way, it becomes a segment of the fourth (or third) 
type. The arrangement ai,bi,a 2 ,b 2 ,... ,an,bn can be converted to another 
arrangement in the following steps: 

1. Extract the uncircled segments of the first and second type from the 
arrangement, modify them as described above, and put the modified seg¬ 
ment back in the same place. 

2. Number the uncircled segments of the third and fourth type from left to 
right. As explained above, they must be equally numerous. 

3. Replace the ith uncircled segment of the third type by the modification 
of the Rh uncircled segment of the fourth type. Similarly, replace the 
zth uncircled segment of the fourth type by the modification of the ith 
uncircled segment of the third type. 

The shuffle vr* is dehned by setting 7r*(i) equal to the ith number in the ar¬ 
rangement constructed in this manner. By construction, vr* is a shuffle from 
(1,2)"' to D 2 . Further, the number of descents of vr* must equal the number 
of descents of vr for the following reason. If there is a descent or an ascent 
between two circled positions in the arrangement ai, 61,02,62 ,..., a„, it 
remains unchanged. Also, the modification of uncircled segments described 
above preserves the number of descents, although it changes their locations. 
Finally, if a circled number is greater than (or less than) a single uncircled 
number, it must be greater than (or less than) all uncircled numbers and, 
therefore, the number of descents between circled and uncircled numbers in 
the arrangement also remains unchanged. 

It is possible to map a shuffle from (1,2)" to D 2 to a shuffle from (1,2)" 
to Di using the same procedure. That map would be the inverse of the map 
dehned above. Therefore, the map from shuffles vr to shuffles -tt* dehned 
above is a bijection. □ 

The converse of the above Theorem 4.1 appears to be true as well. The 
number of equivalence classes seems to be given by the simple formula 
(n-1-3)2"“^. One can attempt to prove this by hnding unique representatives 
for equivalence classes and then counting them as in the proof of Theorem 
3.2. We have derived a method to construct unique representatives for equiv¬ 
alence classes of R. However, we have not yet devised a method to count 
the number of unique representatives. 
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5. Card games. Some inferences about the mixing times for common 
card games such as blackjack and bridge can be drawn using results given 
in the preceding sections. Let 5 be a finite set and let p be a probability 
distribution on S. Then the total variation distance of p from the uniform 
distribution is given by b(®) “ f^l- ^ deck of 52 distinct cards, 

the total variation distance remains close to 1 until the number of riffle 
shuffles exceeds 4. The total variation distance falls below 0.5 when the 
number of shuffles is 7 and this can be taken to be the mixing time [2]. 
Another distance defined in [1] is the separation distance. The separation 
distance of p from the uniform distribution is maxs£5(l — |S'|p(s)). Like the 
total variation distance, the separation distance has a maximum of 1 and 
a minimum of 0. However, it leads to a more demanding notion of mixing 
as the number of riffle shuffles of a deck of 52 distinct cards needed to 
make the separation distance no more than 1/2 is 11. The use of entropy to 
understand mixing is discussed in [14]. The validity and limitations of the 
maximum entropy model of riffle shuffles are discussed in [3] and [5]. 

In the game of bridge, 52 distinct cards are dealt to four players. To apply 
the results proved in Section 2, we need to assume that the hrst 13 cards 
are dealt to one player, the next 13 to another and so on. Let the deck D be 
a permutation of the multiset {1^^, 2^^,3^^,4^^}. Let pD be the transition 
probability from D to 1^^, 2^^, 3^^,4^^ under an o-shuffle. This transition 
probability can be obtained using (2.6). The probability that the first player 
is dealt cards originally in the positions occupied by cards labeled 1 in D, 
that the second player is dealt cards originally in the positions occupied by 
cards labeled 2 in D, and so on after an a-shuffle is equal to pD- Therefore, 
the distance of the probability distribution p^ over decks from the uniform 
distribution will indicate the closeness of a deal after an a-shuffle to a random 
deal to four players. If the separation distance is used to define mixing, an 
application of (2.6) with D = (4,3,2,1)^^ shows that the separation distance 
is greater than 0.5 after 10 riffle shuffles. Therefore, the mixing time is 
still 11 riffle shuffles. The total variation distance involves a sum with a 
great number of terms and the results of Sections 3 and 4 indicate that the 
recognition of equalities among the transition probabilities pz) is unlikely to 
make this sum tractable. However, a Monte Carlo procedure for evaluating 
this sum, which will be described elsewhere, implies that the mixing time 
is 6 riffle shuffles when the total variation distance is used. If the cards are 
dealt to the players in cyclic order, which is the common practice, the mixing 
times will almost certainly be lower. 

In the game of blackjack, the distinction between the suits is ignored. 
We assume the starting deck to be 1^, 2'^,..., 13“^. Application of Theorem 
2.1 and (2.1) shows that the separation distance from the uniform distribu¬ 
tion over decks becomes less than 0.5 after 9 riffle shuffles. Again, a Monte- 
Carlo procedure has to be employed to find the total variation distance. It 
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then follows that the total variation distance becomes less than 0.5 after 
only 4 riffle shuffles. 
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